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Experimental investigation of double layers in expanding plasmas 
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Double layers (DLs) have been observed in a plasma reactor composed of a source chamber 
attached to a larger expanding chamber. Positive ion beams generated across the DL were charac¬ 
terized in the low plasma potential region using retarding field energy analyzers. In electropositive 
gases, DLs were formed at very low pressures (between 0.1 and 1 mTorr) with the plasma expansion 
forced by a strongly diverging magnetic field. The DL remains static, robust to changes in boundary 
conditions, and its position is related to the magnetic field lines. The voltage drop across the DL 
increases with decreasing pressure, i.e. with increasing electron temperature (around 20 V at 0.17 
mTorr). DLs were also observed in electronegative gases without a magnetic field over a greater 
range of pressure (0.5 to 10 mTorr). The actual profile of the electronegative DL is very sensitive 
to external parameters and intrusive elements, and they propagate at high negative ion fraction. 
Electrostatic probes measurements and laser induced photodetachment show discontinuities in all 
plasma parameters (electron density, electron temperature, negative ion fraction) at the DL position. 
The voltage drop across the electronegative DL is about 8 V, is independent of the gas pressure and 
therefore of the electron temperature. 

PACS numbers: 


I. INTRODUCTION 


An electric double-layer (DL) is a region within a 
plasma in which localized charged regions can sustain 
a potential difference and cause electron and ion acceler¬ 
ation [l[. A DL can be thought of as an internal sheath 
separating plasmas at different potentials. Following gen¬ 
eral terminology, in this paper the high plasma potential 
region is referred to as the ’’upstream” region and the 
low plasma potential region across the DL as the ” down¬ 
stream” region. DLs occur naturally in a variety of space 
plasma environments and are of great interest in astro¬ 
physics. In the laboratory, they have been experimen¬ 
tally studied for decades. In most of the experiments, 
the DLs were created in double plasma devices where 
two plasmas were externally biased at independent po¬ 
tentials. Due to the electrical bias, a current is flow¬ 
ing through these current-driven DLs. Since the work 
of Charles and Boswell 0, there has been considerable 
interest in static current free double layers that are self- 
consistently created in low-pressure expanding plasmas. 
These DLs have applications in space propulsion [1, 
and in modeling of astrophysical processes In these 
experiments, a plasma expansion was forced by a strong 
divergence of a high magnitude magnetic field at very 
low pressure. Detailed investigations of these current- 
free DLs have been performed using electrical measure¬ 
ments [3, [fil-fTljl. laser induced fluorescence [l3 - [l^ and 
numerical simulations [l3 . [l8l| . An upstream ionization 
instability has recently been observed in one of the exper¬ 
imental systems [T^. Despite the strong interest shown 
in this area, no general theory has been developed. A 
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recent model of DL formation has been constructed by 
Lieberman and co-workers and compared with experi¬ 
ment [^ . Chen showed that, based on a one dimensional 
model, the physical process is the same as for single layer 
formation, however no general explanation of how the 
single layer develops into a DL is given [2lj |. 

In this article, we present an experimental study of 
the formation of a DL in a low pressure radio fre¬ 
quency driven electropositive plasma where its expansion 
is forced by a strongly diverging magnetic field. The con¬ 
ditions for DL formation and the evolution of the plasma 
parameters are analyzed using electrical measurements 
and compared to other investigations. Throughout the 
text these DLs are termed ’’electropositive DLs”. 

In a recent letter, we showed that a DL could also 
be created in an expanding plasma due to the presence 
of negative ions, without a static magnetic field, and 
at higher pressures [131. Electronegative plasmas are 
widely studied not only because of their relevance for 
processes in the microelectronic industry, but also be¬ 
cause their fundamental properties are significantly dif¬ 
ferent from electropositive discharge plasmas [lH chap. 
8]. Some of the most important and interesting prob¬ 
lems are related to the existence of various kinds of non¬ 
linear potential structures and instabilities. Nonlinear 
potential structures occur in plasma stratification be¬ 
tween an electropositive outer shell (no negative ions in 
the plasma) and an electronegative core (negative ion 
containing plasma). This feature was theoretically stud¬ 
ied and experimentally observed [13, 1131 • Fur¬ 

thermore, it has been shown that inductive discharges 
are subject to relaxation oscillations near the capacitive 
to inductive mode transition intermediate powers) in 
the presence of negative ions (3ll - l34| . At sufficiently high 
power (in the inductive mode), the plasma is not subject 
to these relaxation oscillations. However, in an expand- 
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ing plasma, a periodic propagating phenomena had been 
identified by Tuszewski and co-workers [H, . Plihon 

and co-workers [s^ showed that this periodic propagat¬ 
ing structure is a double layer for their experimental con¬ 
ditions. We remind here the basic features of this phe¬ 
nomenon in Ar/SFe gas mixtures. When the SFe con¬ 
centration is above 13%, the discharge is non stationary, 
with plasma parameters fluctuating due to the periodic 
formation and propagation of a double layer occuring at 
the interface between the source and the diffusion cham¬ 
ber. For lower SFs concentrations (namely between 8 
and 13% at 1.5 mTorr), the double layer appears to be 
static, and eventually disappears for SFg concentrations 
below 8%. The present article details an extensive ex¬ 
perimental study of the static electronegative DL, in the 
same experimental device as the one used for the study 
of electropositive, magnetized DLs. Ar/SFe and Ar/02 
gas mixtures are investigated. Although the physical pro¬ 
cesses leading to the formation of the electronegative DLs 
are different from those involved in the formation of the 
electropositive DL, we provide a comparison for the two 
cases. 

This paper is divided into five parts organized as fol¬ 
lows. Part II presents the experimental setup and diag¬ 
nostics used to determine the plasma parameters. Part 
III presents a parametric study of the electropositive DL 
in a low pressure (below I mTorr), strongly magnetized 
plasma, and emphasizes similarities with previously pub¬ 
lished work. Part IV gives a detailed experimental inves¬ 
tigation of the electronegative DL in Ar/SFg and Ar/02 
gas mixtures. Some conclusions are given in part V. 


II. EXPERIMENTAL SETUP AND 
DIAGNOSTICS 


A. The plasma source 

The plasma source is shown schematically in figure [T] 
It consists of a source chamber on top of a 32 cm diam¬ 
eter diffusion chamber. The source is a 15 cm diameter, 
30 cm long and 0.9 cm thick pyrex cylinder surrounded 
by a double saddle field type helicon antenna [s^. The 
fan-cooled antenna is powered through a close-coupled 
L-type matching network by an rf power supply oper¬ 
ating at 13.56 MHz and capable of delivering up to 2 
kW forward power. The time-averaged input power was 
recorded as the difference between the time-averaged for¬ 
ward and reflected powers. The pyrex cylinder is housed 
in a 20 cm diameter and 30 cm long aluminum cylinder. 
A grid attached to the other end of the source tube con¬ 
fines the plasma and isolates it from a turbo-molecular 
pump that routinely maintains base pressures of 10“® 
mbar. The boundary conditions terminating the plasma 
can be modified, as specified in section IIIIDI The dis¬ 
charge was operated in pure Ar, Ar/SFe and Ar/02 gas 
mixtures. The partial gas pressures of Ar and the molec¬ 
ular gases were determined by controlling the flows and 


the total pressure was measured using a Baratron pres¬ 
sure gauge. 

A static magnetic field can be generated by a d.c. cur¬ 
rent in two coils surroundinfg the source tube. The coils 
are designed such that the magnetic field is axial and uni¬ 
form in the source and strongly divergent at the plasma 
source - diffusion chamber junction. Figure [T] shows the 
normalized axial magnetic field magnitude (solid line) 
and the derivative (dots). Results presented in section Hill 
were obtained with a d.c. magnetic field in a pure ar¬ 
gon plasma, while in section HVl the reactor was operated 
without a d.c. magnetic field. The magnitude of the 
magnetic field mentioned in the text is the maximum 
field obtained along the source axis. 


B. Langmuir probe operation 


The measurements reported here were performed along 
the revolution axis {z axis) of the discharge. All probes 
can also be inserted in the mid-plane of the diffusion 
chamber and radial dependence of the plasma parame¬ 
ters can be measured. The plasma parameters were de¬ 
termined using a passively compensated Langmuir probe 
[ 3 ^ . of 0.25 mm diameter and 6 mm long platinum wire 
tip. The plasma potential, electron density and electron 
temperature were deduced from the I(V) characteristics 
of the cylindrical probe using a Smartsoft data acquisi¬ 
tion system (40j | . The plasma parameters were also deter¬ 
mined from usual I-V curve processing (^ : the plasma 
potential from the zero of the second derivative of I-V 
traces, the plasma density from the electron current at 
the plasma potential, and the electron temperature from 
a linear fit of the logarithm of the electron current (as¬ 
suming a maxwellian electron distribution function). The 
electron energy probability function was determined from 
the second derivative of the I-V curves [ 2 ^, p. 189] 
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where e is the electric charge, me the electron mass, 
Ap the probe collection area and Ig the electron cur¬ 
rent. The electron current was taken to be the difference 
between total collected current and positive ion satura¬ 
tion current, which is assumed to scale as the square 
root of the biasing voltage relative to the plasma po¬ 
tential [12 . According to [4l| we computed the sec¬ 
ond derivative of the I(V) curves using a Savitzky-Golay 
(S.G.) filter [ 4 ^ of order 3, over 7 points. Since the filter 
works only for equally spaced voltage points, we com¬ 
puted the derivatives with respect to the supply voltage 

and used — ("l — R-rrr ^—) H to cal- 

■iVprobe ‘^Vupply V ^Vupply / 

culate the derivatives with respect to the probe voltage. 
The electron density was also obtained from the integra¬ 
tion of the eepf 
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An effective electron temperature was determined from 
the mean electron energy calculated from the eepf 

2 /*+00 

Tes = T. - / eepi{V)V^^^dV (3) 

When operating with the molecular gases, the plasma 
can contain a large number of negative ions. It is possi¬ 
ble to determine the negative ion fraction a, defined as 
the negative ion density over the electron density ratio, 
using a two probe technique . This technique is based 
on the comparison of the Langmuir probe current at the 
plasma potential (the electron thermal current) with the 
positive ion saturation current from a planar probe, ob¬ 
tained in both an electropositive and an electronegative 
plasma. The calculation relies on a sheath model devel¬ 
oped in [23 and requires an estimate of the negative ion 
temperature. 

C. Retarding field energy analyzer 


ion distribution function is centered around the plasma 
potential with a dispersion due to the ion temperature or 
the RFEA resolution (in our case the RFEA resolution 
is such that the ion temperature cannot be measured). 
Eor measurements downstream of a DL, a double peaked 
lEDE is expected, with one peak centered at the plasma 
potential, and the other corresponding to an ion beam; 
the positive ion saturation current is thus the sum of (i) 
an ion beam current from ions created upstream and ac¬ 
celerated within the DL, and (ii) a local ion current from 
ions created in the downstream plasma. In order to probe 
the temporal development of the lEDF when pulsing the 
plasma, we developed a time resolved data acquisition 
system. The discriminator grid biasing voltage was kept 
constant during 256 plasma pulses and increased step by 
step afterward. The collected current was recorded us¬ 
ing a large memory digital oscilloscope, averaging over 
256 pulses to reduce noise (the temporal behavior was 
checked to be reproducible between each pulse). Voltage- 
current characteristics were then reconstructed and the 
time resolved lEDFs were computed. 


In order to investigate the positive ion distribution 
function, we have developed a retarding field energy an¬ 
alyzer (RFEA) consisting of four grids and a collector 
plate |4fil - l4^ . The four grids are made of nickel wires of 
11 microns in diameter spaced by 40 microns; each grid 
had a 60 percent transparency. The entrance grid was 
grounded, the second grid was biased at -50 V to repel 
the electrons, the third grid was used to select the ion 
energy by scanning the voltage from -10 to 60 V, and 
a fourth grid was inserted before the collector in order 
to minimize the effect of secondary electrons. The grids 
were spaced by 0.25 mm so that the total system length 
was about one millimeter. When assembled, the ana¬ 
lyzer was 35 mm long by 50 mm in diameter and the 
plasma particles enter the analyzer through a 2 mm hole 
in a 0.3 mm thick stainless steel orifice plate. The plate 
is in electrical contact with the analyzer housing, which 
is connected to the grounded diffusion chamber of the 
reactor. The RFEA is supported by a 3 cm diameter 
aluminum tube that is differentially pumped and can be 
moved along the vertical axis. The ion current to the 
collector Ic was recorded as a function of the dc voltage 
applied to the discriminator Vd and is given by: 


IM = cArG^ 
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( 4 ) 


where Ar is the collecting surface, G the grid trans¬ 
parency, Vo the speed of the ion accelerated by the dis- 
crimator grid at potential Vd, fiv) the ion distribution 
function and M_|_ the positive ion mass. The power 4 of G 
accounts for the transmission across the four grids. The 
ion energy distribution function (lEDF) was computed 
as When operating in a low pressure plasma, the 


D. Probe based laser induced photodetachment 


When operating the discharge in electronegative gas 
mixtures, the negative ion fraction a was measure d by 
probe based laser induced photodetachment (LIP) 

The experimental arrangement consisted of a pulsed 
cylindrical laser beam aligned to be collinear with an un¬ 
compensated probe positively biased. The negative ion 
fraction is determined from measurements of the electron 
saturation current in the absence of the laser pulse (Ig dc, 
since the laser duty cycle is less than 10“^) and the in¬ 
crease in the current (A/e) immediately after the pulse, 
using the relationship 


A/e n_ 

d^e dc Vlg 
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where n_ and rie are the negative ion an electron densi¬ 
ties, respectively. The necessary conditions for the 
laser beam diameter (set to be 6 mm), the laser beam 
power and probe voltage (set to be -|- 35 V) are set such 
that the negative ion fractions were measured correctly. 
The laser photon energy required to detach an electron 
depend on the electron affinities of the molecular gas used 
(SFe or O 2 ). 

When operating in Ar /02 mixtures, 0“, and O;^ 
can be produced in high density low pressure plasmas, 
with electron affinities of 1.46 eV, 0.44 eV and 2.1 eV, 
respectively [10, HIl. Thus the photon energy from a 
frequency doubled Nd:YaG laser (532 nm - 2.33 eV) is 
sufficient to photodetach all negative ions. 

When operating in Ar/SFe mixtures, negative ions are 
produced from collisions involving SFe fragments, and 
F“, SF /, SF^ and SFg are the dominant negative ion 
species [52-1^ . The electron affinity for these species are 
3.4 eV O, 2.88 eV [1^, 3.79 eV [13 and 3.16 eV [13 
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respectively. However the electron affinity of SF 5 is still 
debated and might be as much as 4.2 eV [^. The use of 
a frequency-tripled Nd:YaG (355 nm) provides photons 
with sufficient energy to detach most of the negative ions 
(3.5 eV). It is possible that SF^ is not photodetached and 
this will cause discrepencies with negative ion fraction 
determined by the two probe technique. 

The laser power dependence of the LIP signal is given 

by 

A/e , ( aPAt\ 

where hv is the photon energy, P is the laser beam power, 
At the laser beam pulse length, is the laser beam sec¬ 
tion, and (T the LIP cross section for photodetachment of 
a negative ion at a particular laser energy (a is zero for 
energies hv lower than the electron affinity). From pub¬ 
lished cross sections, it is possible to identify the dom¬ 
inant negative ion in the discharge by fitting the beam 
power dependence of the LIP signal to theoretical depen¬ 
dence. 


III. DOUBLE LAYERS IN LOW PRESSURE 
MAGNETIZED ELECTROPOSITIVE PLASMAS 

In this section results obtained in an Ar plasma, at 
very low pressure and with d.c. current flowing in the 
source coils are presented. The measurements confirm 
the existence of DLs in these types of discharges and 
complements previous work. 


A. Existence of a double layer and spatial 
dependencies 

The RFEA was inserted into the diffusion chamber, 
facing the source. A typical current-voltage characteris¬ 
tic obtained from the RFEA is given in figure [5] (solid 
line) when operating at sufficiently low pressure (below 
1 mTorr) and moderate magnetic field (above 60 G). In 
order to provide a systematic analysis, the experimental 
lEDF (dashed line in figure [2|) is fitted as the sum of 
two independent gaussian functions according to a least 
square procedure (dotted line). 

The gaussian functions are centered around Vp and 14 
and their integrals are referred to as Ip and /{,. The 
determination of the plasma parameters is processed as 
follows: (i) the plasma potential is the position of the 
first maximum I 4 , (h) the local ion flux (Fp) is the to¬ 
tal flux (F+) times the ratio Ip/{Ip + h), (hi) the beam 
potential is the position of the second maximum I 4 , (iv) 
the ion beam flux (F;,) is the total flux (r+) times the 
ratio h!{Ip + h)- The total flux r+ is deduced from the 
collected current when the discriminator is set to 0 V, 
corrected from the grid transparency according to equa- 
tion[31 The total flux measured by the RFEA was verified 
against planar probe flux measurements. The potential 


drop across the DL is computed a.s ^dl = Vh— Vp. The 
first peak of the lEDF corresponds to ions that are ac¬ 
celerated at the Bohm velocity VBohm when entering the 
sheath in front of the RFEA; the local plasma density is 
computed as 


cG^UBohm 


( 6 ) 


The second peak of the lEDF is related to ions created 
upstream of the DL that are accelerated within the DL. 
The beam velocity is computed as 


Vb 


l 2q{Vb-Vp) 

M+ 


and the beam density as 
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Figure |3] shows results obtained at 0.17 mTorr, 90 G, 
250 W as a function of axial position. For positions be¬ 
low 28 cm, the lEDF clearly shows two peaks when the 
RFEA is facing the source while one single peak is present 
above that axial position. When the RFEA is rotated by 
90°, only one peak is observed, centered at the same po¬ 
tential as the first peak when the RFEA is facing the 
source. This indicates the presence of an abrupt increase 
of the plasma potential at z Ri 28 cm, associated with 
a dramatic increase in the plasma density, i.e. the dou¬ 
ble layer. This DL is associated with the presence of an 
ion beam in the low potential region, the speed of which 
is determined by the potential drop across the DL. The 
plasma potential and plasma density are roughly con¬ 
stant below 28 cm, with the presence of a beam at 42 
V, the amplitude of which decreases exponentially when 
moving away from the DL; above 28 cm, one single peak 
at about 42 V is obtained, with a local ion density 10 
times higher than downstream (excluding the ion beam 
density), and increasing when moving up into the source 
tube. The plasma density is determined using an elec¬ 
tron temperature assumed to be spatially homogeneous, 
which was measured to be 5±0.5 eV at 0.17 mTorr in 
the diffusion chamber. With the voltage difference be¬ 
tween the beam and local plasma component being 13 
eV, the beam velocity is -v/SuBohm ~ 2.2wBohm according 
to equation [71 The beam measured downstream of the 
DL is therefore supersonic. 

It should be noted that only indirect measurement of 
the DL was performed (via interpretation of the lEDFs 
and beam measurements). It was not possible to measure 
the plasma potential discontinuity by using the Langmuir 
probe. The probe disturbed the plasma such that the 
matching had to be dramatically changed when crossing 
the DL (the mismatch could reach 7 % of the coupled 
power when moving the probe by only 1 cm a few cen¬ 
timeters when positioned downstream of the DL). 

These features are very similar to the work presented 
in references The most significant difference 
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from the above references is the spatial variation of the 
ion beam flux in the downstream region. The main phys¬ 
ical process involved in the damping of the ion beam is 
ion-neutral charge exchange collisions. The charge ex¬ 
change mean free path for thermal ions in low pressure 
plasmas (the ion temperature being around 0.1 eV) is 
[cm] = 3/p where p is the pressure in mTorr [^ . 
The velocity dependence of the charge exchange cross 
section j2^ leads to introducing a correction factor for 
the mean free path for a 15 eV beam and to estimating 
the effective mean free path as A = Ai/0.7 [1^ (at 0.17 
mTorr, A = 25cm). Although the actual geometry is a 
complex cylindrical geometry, we assume a one dimen¬ 
sional descripion of the ion beam density for the sake of 
simplicity. In order to compare the experimental beam 
decay to models, assumptions on the actual DL 3-D pro¬ 
file ar required since direct measurement of the profile 
is not easily achieved. We therefore consider two cases, 
namely a flat profile (perpendicular to the z axis) fol¬ 
lowing the statements given in Q and a hemi-spherical 
profile. These profiles are sketched in figure SKa). As¬ 
suming a flat profile, the one dimensional ion beam flux 
(T+b) conservation is 

rVb = -r+b/A 

where the prime denotes the first derivative with respect 
to z. Integration of this equation with Zq = 28 cm 
leads to an ion beam density spatial evolution n^b(z) = 
n+boexp(—(Z q — z)/\), plotted as a dashed line in fig¬ 
ure [3Kb). The agreement with the experimental values is 
not satisfactory. The uncertainty over the exact pressure 
value, discussed in the following, could lead to a lower 
actual mean free path. Another reason may be due to 
the actual 3-D profile of the DL. Assuming the DL has a 
spherical shape, the ion beam flux conservation in spher¬ 
ical coordinates is 


r-\r^T+b)' = -r+b/A 


where the prime denotes the first derivative with respect 
to r (the origin of the spherical coordinates being at point 
S in figure HKa)). Integration of this equation from the 
DL position assuming a Rq radius for the DL profile leads 
to an ion beam density spatial evolution 


f r- Ro 

n+b(r) = n+bo-^exp I-^— 


(9) 


which can be expressed as a function of z 


n+b(z) = n+bo 


Rl 
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exp 


Zo - z 
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( 10 ) 


This theoretical evolution is plotted on figure ISKb) as a 
solid line, and agrees reasonably well with the experimen¬ 
tal values, especially above 20 cm. 

Understanding the spatial damping of the ion beam is 
important if this device is to be used as a space thruster 
as was proposed earlier Further investigations on 


the 3-D profile of the double layer is needed and opti¬ 
cal diagnostics could provide useful informations if the 
experiment could be reproduced in a device with full op¬ 
tical access to the source tube. Another important is¬ 
sue is the actual value of the charge exchange mean free 
path.; precise knowledge of the operating pressure is thus 
needed. The operating pressure measured is the outer re¬ 
gion of the diffusion chamber, using a 100 mTorr range 
Baratron pressure gauge. The relative error in the range 
0.1-1 mTorr is the Baratron measurement could be as 
high as 50%, even though absolute pressure values for 
DL existence are similar to the one obtained by Charles 
in 


B. Magnetic field and pressure dependence 

Similarly to previous studies led at ANU [1,[ll| , and 
at WVU and Princeton University [isl - fT^ . two impor¬ 
tant experimental control parameters were identified in 
the DL formation process: the magnetic field magnitude 
and the operating pressure. This subsection presents a 
study of the DL sensitivity to these parameters. 

The first observation is that a sufficiently high and di¬ 
verging static magnetic field is required to form a DL. 
Figure jSKa) shows the behavior of the DL as a function 
of the magnetic field strength at a pressure of 0.17 mTorr, 
and a power of 250 W. For magnetic field magnitudes be¬ 
low 45 G, the lEDF shows only one peak, with a fairly 
high plasma density (10^®m“^). Above 45 G, the dou¬ 
ble layer is formed. While the plasma density decreases 
slightly with increasing magnetic field (from 3 10^®m“^ 
at 45 G to 1 10^^m“^ at 90 G), the plasma potential 
increases slightly from 28 V at 45 G to 30 V at 180 G. 
The potential drop across the DL (around 15 Af for these 
conditions) varies weakly with the magnetic field mag¬ 
nitude; these variations remain however within the error 
bars of both the current measurement and the systematic 
lEDF processing. The beam density variations appears 
to follow the local plasma density behavior. As in previ¬ 
ous studies, the DL characteristics are fairly independent 
of the magnetic field above a critical value. However the 
actual shape and amplitude of the DL seems to be re¬ 
lated to the actual geometry of the system and magnetic 
field topology. Topological effects might explain the dis¬ 
crepancy between the high potential drop across the DL 
observed by Gharles 0 and the lower potential drop ob¬ 
tained in this work and in Ref. [l^. The DL position 
also depends on the geometry of the reactor, but remains 
close to the source - diffusion chamber junction (slightly 
downstream in WOMBAT , and slightly upstream in 
Ghi-Kung 0, HELIX [T^ and our source). These ex¬ 
periments show that above the critical magnetic field at 
which the DL forms, the DL characteristics do not de¬ 
pend on the magnetic field magnitude; thus the topology 
of the magnetic field might be the key parameter. The 
position of the DL in all these systems is close to the 
position of the maximum in the magnetic field gradient. 
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The second observation is that the double layer only 
appears at sufficiently low operating pressures. Fig¬ 
ure Eljb) shows the influence of the pressure on the DL 
characteristics for a 90 G magnetic field and an operating 
power of 250 W with the RFEA positioned at z = 20 cm. 
No discharge could be sustained below 0.08 mTorr and 
between 0.08 and 0.1 mTorr, the measured lEDFs show 
only one peak, i.e. no DL is present. Between 0.1 and 
1 mTorr, two peaks are clearly observed, providing evi¬ 
dence of a DL. When increasing the pressure, the plasma 
potential downstream of the DL remains constant, while 
the upstream potential decreases. The DL potential drop 
decreases dramatically with increasing pressure, as was 
determined in the model by Lieberman et al. (20|. The 
beam flux also decreases with increasing pressure due to 
the decrease in acceleration from the DL potential drop 
and the increase in collision frequency. The plasma den¬ 
sity scales as expected at constant power: increasing with 
increasing pressure, since the electron temperature re¬ 
quired to sustain the discharge decreases. At pressures 
above 1 mTorr the two peaks start to overlap. The reso¬ 
lution of the RFEA does not allow to distinguish between 
an upper pressure threshold for DL formation or a grad¬ 
ual decay of the DL potential drop to zero. In any cases, 
the DL disappeared at 3 mTorr. 

The DL dependence on pressure and magnetic field 
strength are similar to previous experimental investiga¬ 
tions and are correctly caught by the models developed 
in and [2l| . Chen based his model on the hypothesis 
of a plasma frozen to the magnetic field lines [^. The 
magnetic confinement can be estimated from the value of 
the ratio Wce/^'m [13 j where ujce is the electron cyclotron 
pulsation (proportional to the magnetic field magnitude) 
and Vm is the total electron collision frequency (propor¬ 
tional to the pressure). Thus the assumption of a plasma 
frozen to the field lines is satisfied at high magnetic field 
and/or at low pressure, which are the conditions for DL 
occurrence. In this model, the potential drop across the 
DL is dependent on the electron temperature and the 
condition for development of a non-neutral region is for 
the plasma column radius to increase by 28%. An anal¬ 
ysis of the calculated magnetic field lines (assuming a 
plasma frozen to the field lines) leads to a DL position 
at z Ri 30 cm according to this condition, which is close 
to the observed value (28 cm) and to the maximum of 
the magnetic field gradient. The model by Lieberman et 
al. [13 explicitly accounts for the pressure dependence of 
the DL. All trends emphasized in this study are caught 
by the model, namely the pressure dependence of the po¬ 
tential drop across the DL , a low pressure limit for DL 
existence due to a low ionization term from accelerated 
electrons upstream, and a high pressure limit when no 
ionization is needed from accelerated electrons upstream. 
However, the model in [13 does not include a magnetic 
field, which is a key parameter of this experiment. 

The experimental investigations also show that the 
potential drop across the DL was fairly insensitive to 
the coupled power. However the ion plasma density 


to ion beam density ratio in the diffusion chamber is 
highly dependent on the power, increasing with increas¬ 
ing power [^ . 

C. Ignition of the electropositive double layer 

Using the time resolved data acquisition system, the 
ignition and the time development of the DL was mea¬ 
sured. The discharge was operated in a pulsed mode, 
with a 750 /xs on period and 500 fis off period. Although 
the lifetime of some Ar metastables states is of the or¬ 
der of hundreds of microseconds at the current operating 
pressures, the ’off’ period allows sufficient time to destroy 
all charged species. Time-resolved lEDFs are shown in 
figure UKa) but need to be analyzed bearing in mind the 
generator time scale. Analysis of the pulsed mode oper¬ 
ation of the radiofrequency power generator on a 50 H 
load showed that the typical generator rise-time is 150 
fis. The transient plasma operation time will thus last at 
least 150 /xs. 

The lEDFs are processed to follow the amplitude and 
potential of the two peaks, which are given in figurelBKb). 
The beginning of the ’on’ time is located at t = 20 /xs, 
and only one peak is present before t = 50^s. The beam 
potential reaches its steady-state value over 10-20 /xs, 
while the downstream plasma potential keeps on decreas¬ 
ing (from 38 to 30 V) over the next 50 p,s. Formation of 
the double layer occurs less than 30 /xs after plasma break 
down and the steady state potential values are reached 
well before steady state density is reached. The ampli¬ 
tude of both the beam component and local plasma com¬ 
ponent increase during the first 150 /xs of operation, and 
increasing at the same rate. 

These observations confirm and complement Charles’s 
investigations [6l|, by providing a detailed analysis of 
the lEDF temporal evolution. The time scale for the DL 
to develop is found in these experiments to be around 
100 /xs. However,in these experiments, an exact evalua¬ 
tion of the DL time development was not possible since 
the ignition time is less than the generator risetime (the 
overshoot in figure 4 of reference [^ might be due to 
the generator temporal dynamics). Temporal evolution 
of the downstream beam velocity from LIF in the HELIX 
system [l3 over long pulses showed that a steady state 
was not reached after 100 ms of operation. In our sys¬ 
tem, it may be that the DL forms and develops rapidly 
over the first 200 /xs, and then experiences a slower rate 
of change as was observed in the HELIX reactor. In our 
system, it was not possible to operate with longer pulses 
so as to observe a change in the temporal dynamics of 
the DL on large time scales (hundreds of ms to seconds). 

D. Influence of the boundary conditions 

We have investigated the influence of changing the 
boundary conditions at the end of the source tube. Three 
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different configurations were experimentally studied: (a) 
an insulating grid positioned 2 cm above the top of the 
pyrex source tube with a dc grounded metal cylinder in- 
between, (b) an insulating grid 2 cm above the top of 
the pyrex source tube with an insulating pyrex cylinder 
in-between and (c) a non-dc grounded conducting grid at 
the top of the pyrex source tube. These configurations 
are sketched in figure [71 All results presented in the pre¬ 
ceding subsections and in the following section were ac¬ 
quired with condition (a), namely an insulating grid at 
the end of the pyrex tube with a piece of grounded metal. 

The bottom part of figure [7] shows lEDFs normalized 
to one at the plasma potential, when the RFEA is po¬ 
sitioned at z = 20cm, for a 0.17 mTorr, 90 G, 250 W 
plasma. The variation of the ion saturation flux is within 
3% for the three configurations presented. The beam am¬ 
plitudes decrease by a few percent (8% from geometry (a) 
to geometry (c)). A noticeable variation of the potentials 
is observed: the local plasma potential decreases: 26.5 V, 

25.5 V and 24.5 V, while the beam potential increases: 

42.5 V, 44.2 V, 46.5 V, for geometry (a), (b) and (c) 
respectively. The broadening of the low potential peaks 
increases slightly when scanning from geometry (a) to ge¬ 
ometry (c). Due to the low operating pressure and large 
mean free path, the RFEA collecting length is large (ap¬ 
proximately one mean free path). The lEDF’s first peak 
broadening could originate as a result of smaller plasma 
potential gradient lengths in the downstream region from 
geometry (a) to (c). A gradient length less than the mean 
free path leads to a non-mono-energetic beam and peak 
broadening. The potential drop across the DL changed 
noticeably with the boundary conditions, the higher po¬ 
tential drop across the DL being obtained for the floating 
conducting grid. This investigation confirms the influ¬ 
ence of the geometry, but these results are significantly 
different from those published in Ref. Q which showed 
that the beam component relative to the plasma compo¬ 
nent could triple when changing the position of the end 
plate by only 2 cm. 

Interestingly, the double layer was formed for all 
boundary conditions under the same experimental con¬ 
ditions, including the partially dc grounded case (con¬ 
dition a). DLs with conditions (b) and (c) have to be 
current free, while DL with condition (a) could be cur¬ 
rent carrying. Whatever its nature (current-free or possi¬ 
bly current-carrying), the DL’s characteristics are weakly 
modified as emphasized in Ref. [2^. Experiments show 
that the same trends are maintained as a function of pres¬ 
sure, magnetic field magnitude and position whatever the 
source boundary conditions. 


IV. DOUBLE LAYERS IN ELECTRONEGATIVE 
MIXTURES 

Double layers have been observed in the same experi¬ 
mental device operating without a magnetic field in elec¬ 
tronegative gas mixtures [12, H^l- In this section, we 


present a detailed experimental investigation of the static 
electronegative DL introduced in [121. Results were ob¬ 
tained in a stable inductive mode, i.e. far from the relax¬ 
ation oscillations between the capacitive and inductive 
modes. The section is organized as follows: we first de¬ 
velop results obtained in Ar/SFg gas mixtures and then 
present results obtained in A 1 /O 2 gas mixtures. 


A. Conditions for DL formation 

We emphasized in the introduction that DL formation 
was highly dependent on the SFg concentration: a DL 
forms at the interface between the source tube and the 
diffusion chamber when the SFg concentration is above 
a critical concentration. This DL remains static for a 
narrow range of SFg concentration (namely 8 to 13% at 

1.5 mTorr), before a propagating mode is established. 
Variation of the SFg concentration within the [8-13%] 
window has a dramatic impact on the DL characteristics 
and overall plasma equilibrium. At a typical 1.5 mTorr 
plasma pressure, the position of the DL on axis is mea¬ 
sured at z = 25 cm for a 8% SFg concentration and as 
low as 18 cm for a 11% SFg concentration. Direct mea¬ 
surement of the 3-D profile is not easily achieved. How¬ 
ever from visual observation through an optical window 
in the diffusion chamber, the DL profile can be seen as 
an interface between a bright plasma in the source and 
a more diffuse plasma in the diffusion chamber. The 
actual 3-D profile has a hemi-spherical shape entering 
the diffusion chamber and is attached to the bottom of 
the source tube. This is similar to the plasma-sac pre¬ 
viously described in Ref. [62] . A change of the profile 
observed when increasing the SFg concentration is given 
schematically in figure |H[b) . The plasma equilibrium is 
also strongly modified with a change in SFg concentra¬ 
tion: the downstream plasma potential is decreasing from 
18 to 15 V for SFg concentrations increasing from 8 to 
11%. The DL’s characteristics are sensitive to any pa¬ 
rameter change; slight changes in pressure or power lead 
to a modification of the DL profile at a given SFg con¬ 
centration. Over time (month to month) similar slight 
changes were observed due to wall contamination or ther¬ 
mal drift for the same SFg concentration. However the 
same trends and values were reproduced over a time pe¬ 
riod of two years. All SFg concentrations specihed in this 
paper are calibrated according to a reference average con¬ 
centration. Furthermore, intrusive probing of the plasma 
parameters can modify the DL position. At low plasma 
density (coupled power around 100 W), the DL formation 
was observed, but its shape could be modified by moving 
the Langmuir probe in its vicinity. This behavior was 
less extreme at higher plasma densities, since the probe 
area is small compared to the plasma section. Inserting 
the RFEA in the diffusion chamber changes the 3-D pro¬ 
file of the DL. When the RFEA is moved upward in the 
chamber, the DL is ’’attracted” towards the RFEA. It 
was not possible to measure the lEDF close to the DL 


without disturbing it: for a 8% SFg concentration, the 
DL profile is not modified (crossing the z axis at 25 cm) 
when the RFEA is at z = 8-9 cm, , while the RFEA 
destroys the DL when positioned at z = 15 cm. 

All results presented in the following subsections were 
obtained slightly above the critical SFg concentration for 
DL formation, i.e. the DL profile is close to the 8% SFe 
case plotted in figured The reader should keep in mind 
the high sensitivity of the DL to intrusive probes, espe¬ 
cially when analyzing RFEA results. 

Electronegative DL formation was observed to be eas¬ 
ier than electropositive DL formation. The electronega¬ 
tive static DL was created for all pressures in the range 
0.5-6 mTorr, for all power from 50 to 1500 W (i.e. both 
in the capacitive and inductive power coupling modes). 
For any power-pressure range, a static electronegative DL 
forms within a narrow window of low SFg concentration. 
At higher pressures (above 6 mTorr), the electronegative 
DL was always observed propagating. 


B. Spatial evolution of electron parameters 

Figure [8] shows the spatial evolution of the plasma pa¬ 
rameters on axis for two SFg concentrations: 4% (dots) 
and 8% (open circles). The lower SFg concentration pro¬ 
files correspond to the no DL case and are similar to pure 
argon prohles (when operated without a d.c. magnetic 
field); the plasma is created in the source and diffuses 
continuously into the diffusion chamber. The profiles for 
a 8% concentration clearly show a discontinuity at the 
interface between the source and the diffusion chamber; 
the DL is created at this interface (seen at z Ri 25cm 
for all plasma parameters). In the source, the plasma 
parameter prohles are similar in both the non-DL and 
DL cases. All results presented in this subsection and 
in subsections IIV Cf IIV PI were obtained for a 1.1 mTorr 
pressure measured when the discharge is ”ofF’. Due to 
dissociation of the molecular gas, the pressure lies in the 
range 1.4-1.5 mTorr when the plasma is ignited. 

The potential drop across the DL is difficult to measure 
from the Langmuir probe axial prohles, since upstream 
of the DL a continuous decrease in potential, similar to 
a pre-sheath prohle, is observed, leading to a somewhat 
indeterminate beginning of the non-neutral region. How¬ 
ever this potential drop can be estimated to be around 7 
V. The potential drop can also be estimated from equa¬ 
tion 6 in Ref. [l^ which estimates the DL amplitude 
from an electron hux balance between the upstream and 
downstream plasmas, and gives a value of 6.7 V (when 
assuming the upstream plasma edge at z = 25 cm). 

Upstream and downstream of the DL, the electron dy¬ 
namics is clearly different. Figure [HDa) shows the eepf as 
a function of energy and axial position. The DL separates 
maxwellian-like distributions downstream and truncated- 
maxwellian like distributions upstream. Figure[SDb) gives 
a view of eepfs at axial position 22 cm and 28 cm (namely 
3 cm downstream and 3 cm upstream of the DL), plot¬ 


ted as symbols. The upstream distribution is similar to 
those obtained in low pressure inductively coupled plas¬ 
mas [6^. The truncature energy Cg tends toward the 
sheath voltage (difference between the plasma potential 
and floating potential in our situation) at very low pres¬ 
sure (below 1 mTorr) and toward the excitation collision 
energy threshold (above 10 mTorr). In figure [HKb) , Cg is 
estimated to be 17 V (slightly higher than the excitation 
collision energy threshold of 11.5 eV in argon). 

The electron density and effective electron temperature 
determined from the eepf are consistent with the results 
presented in hgure [S] The electron density obtained by 
integration of the eepf follows the same trend as the elec¬ 
tron density computed from the I-V curve analysis with 
absolute value discrepancies less than 10%. The spatial 
evolution of the effective electron temperature exhibits 
an abrupt drop at the DL position. However, the effec¬ 
tive electron temperature values differ from those derived 
from direct I-V curves analysis: downstream of the DL 
the effective electron temperature is measured to be 3.8 
- 4 eV, and 5.5 - 6 eV upstream. The abrupt difference 
in effective temperature is in the range 1.5-2 eV. 

The change in the electron dynamics across the DL 
could originate in the shift of the upstream eepf through 
the retarding potential of the DL. A 8.5 V Boltzmann 
shift of the upstream eepf of figure [HKb) is represented 
by the black solid line. This shift cuts off the lower en¬ 
ergy and brings the distribution into agreement with the 
downstream distribution. The 8.5 V value is obtained 
from the difference in potential at z = 22 and z = 28 
cm taken from hgure [5Ka). However the shifted distri¬ 
bution is somewhat different from the downstream dis¬ 
tribution at low energies. This could be due to distinct 
physical processes upstream and downstream (ionization, 
excitation and attachment processes may vary due to a 
difference in the plasma equilibrium and species). The 
electron density and effective electron temperature cal¬ 
culated from the shifted distribution are 8 x 10^^cm“^ 
and 3.9 eV respectively (compared to 9 x 10^®cm“^ and 
3.85 eV from integration of the energy distribution at 22 
cm). 


C. Evidence of a positive ion beam in the 
downstream region 

Similarly to the electropositive case, lEDFs were mea¬ 
sured in the low potential region with the RFEA facing 
the source. Experimental evidence of an ion beam in the 
diffusion chamber, accelerated through the DL, is drawn 
from the lEDFs. Figure [TO] shows an experimental lEDF 
downstream of the DL exhibiting two peaks. Due to the 
RFEA’s resolution, these two peaks are not clearly sepa¬ 
rated. However, as shown in figure [TUI the lEDF can be 
fitted as the sum of two gaussian functions (amplitude 
Ai and deviation ai). Following the analysis presented 
previously for the electropositive DL (section HI), the 
lEDF provides a measure of the plasma potential (1^), 
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the beam speed (Vf, — Vp) and the positive ion flux from 
the local plasma {T +^ 20 - 2 )) and the beam 
flux (r+A 2 C’' 2 /(^icri + j42cr2))- The beam speed is an in¬ 
direct measurement of the potential drop within the DL. 
The dispersion cti of the gaussian centered at the plasma 
potential is smaller than the dispersion 172 of the gaus¬ 
sian centered at the beam potential. The dispersion of 
the peak centered at the plasma potential is due to the 
resolution of the RFEA and is not a measure of the ion 
temperature since CTi is 3 V, which is far higher than the 
typical positive ion temperature in such a plasma (which 
is less than 0.5 eV). Concerning the higher value of <72, it 
is possibme that the ion beam is not mono-energetic and 
spreads out, due to the gradient in the plasma potential 
upstream of the DL (as given in figureUKa)), whose gradi¬ 
ent length is of the same order as the ion mean free path. 
The insert of figure [TT] shows the lEDF acquired with the 
RFEA inserted radially, into the mid-plane of the diffu¬ 
sion chamber at a 6 cm radial extension from the axis. 
This lEDF has only one peak, centered at the plasma 
potential that was determined from the axial lEDF. The 
second peak observed on the axial lEDF is therefore due 
to a population created upstream and accelerated within 
the DL. 

Measuring the lEDF downstream is a non-invasive 
probing method of the DL existence. However, we pre¬ 
viously emphasized that when the RFEA is positioned 
in the upper part of the diffusion chamber, the DL pro¬ 
file was profoundly affected. Keeping this limitation in 
mind, figure [TT] shows the spatial evolution of the plasma 
potential and the beam potential, as well as the local and 
beam fluxes from lEDF measurements. The DL position 
is lowered down to z = 15 cm for the higher probing po¬ 
sitions. Assuming the damping of the ion beam flux to 
originate in ion-neutral charge exchange, the mean free 
path for this process is taken to be, in a similar fash¬ 
ion to the electropositive case, X[cm] = 3/(0.7p[mTorr]) 
where p is the argon partial pressure. The argon pres¬ 
sure was set to be 1 mTorr when the plasma was off. 
Even though the pressure is higher when the plasma is 
struck, we state that this increase in pressure is due to 
SEg molecular dissociation, and that the argon partial 
pressure remains at 1 mTorr, leading to A = 4.3 cm. As 
is similar for the electropositive case, theoretical evolu¬ 
tions of the beam damping are compared to the mea¬ 
sured beam amplitudes in figure [m b') for two assumed 
DL profiles. The flat profile case is plotted as a dashed 
line, while the hemi-spherical profile case (with Rq = 10 
cm) is plotted as a solid line. 

Interestingly, the best fit is obtained for an assumed 
flat profile of the DL, although observations clearly indi¬ 
cate an experimental hemi-spherical profile (with a cur¬ 
vature radius in the range 12-15 cm). The effect of the 
RFEA disturbances were not taken into account in these 
theoretical plots. A decrease of the DL position by the 
RFEA attraction leads to a reduced effective mean free 
path. On the other hand, deeper penetration of the DL 
into the diffusion chamber leads to a reduction of the ion 


flux crossing the DL (since the total flux remains con¬ 
stant) and this effect could compensate for the previous 
one. Damping of the beam is thus weakly understood 
and is theoretically underestimated for the electroposi¬ 
tive DL and overestimated for the electronegative DL. 
Precise measurement of the ion beam damping length 
can only be obtained by the use of non-intrusive Doppler- 
Fizeau shifted laser diagnostics. 


D. Negative ion fraction measurements 


Spatial evolution of the negative ion fraction in the 
diffusion chamber have been presented in reference [ 2 ^. 
The evolution was processed from electrostatic probe 
measurements according to the two probe technique de¬ 
scribed in [ 4 ^. This technique relies on a sheath the¬ 
ory where negative ions are assumed to be in Boltz¬ 
mann equilibrium (the electric force is balanced by the 
thermodynamic pressure which is proportional to the ac¬ 
tual negative ion temperature) as is usually the case in 
low pressure plasmas [64( 1 . It requires the knowledge of (i) 
the dominant positive ion mass, and (ii) the negative ion 
temperature. The negative ion temperature is assumed 
to be proportional to the electron temperature, with a 
coefficient 7 = Te/T_ usually estimated in the range 10 - 
20 in similar discharges [^. Moreover, the resolution 
of the technique is 7 -dependent , but is qualita¬ 

tively limited to a > 2.5 — 3 in similar plasmas. The 
results presented in [22| were processed assuming 7 = 15 
and an effective positive ion mass of 40 a.m.u. These 
results are reproduced in figure IT^ ai . According to this 
technique, the negative ion fraction scales as the square 
root of 7 and the effect of using 7 = 5 is plotted in fig¬ 
ure [T^a) (open circles). The limitations of these results 
should be highlighted before comparison with LIP mea¬ 
surements. The first limitation lies in the negative ion 
temperature estimate, and the validity of a Boltzmann 
relation for negative ions. The second limitation is the 
resolution of the technique: we are not able to measure 
a negative ion fraction lower than 3 for 7 = 15. Hence, 
the actual negative ion fraction value might be in the 
range [0;3] when a value below 3 is measured. Negative 
ion fractions upstream of the DL and in the lower part of 
the diffusion chamber are therefore upper bounds. The 
last limitation is due to the ion beam existence in the 
downstream plasma. The planar probe used is facing 
the source and the positive ion saturation flux has two 
components: a local flux and a beam flux. The sheath 
theory does not apply in the presence of a beam. Since 
electroneutrality reads n+p + n+t, = (I -|- a)ne where n+p 
is the local ion density and n+f, the beam ion density, one 
should use F+p -|- F+f, ^ as the positive ion flux in the 
negative ion calculation, where ub is the modified Bohm 
velocity and Vb the beam velocity. Since Vb ^ ub if the 
negative ion fraction is above 1 , a correction can be esti¬ 
mated by subtracting the beam flux from the total flux. 
Within the first 4 cm downstream of the DL, this leads to 
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a correction factor 0.5 according to the RFEA’s flux mea¬ 
surements shown in figure fTTT b'). We thus believe that a 
few centimeters downstream of the DL the actual value of 
the negative ion fraction from electrostatic probes could 
be half the values plotted in figure fT^ al. 

Probe based laser induced photodetachment (LIP) al¬ 
lows the direct measurement of the negative ion fraction 
as was outlined in section im Figure [T^ bl presents the 
spatial evolution of the negative ion fraction in the diffu¬ 
sion chamber and in the bottom of the source for the two 
SFg concentrations (4 and 8%) previously detailed. Neg¬ 
ative ion fractions measured by the two probe technique 
and displayed in figure fT^ al are for a 8% SFg concentra¬ 
tion (as was emphasized at the beginning of the section, 
the exact position of the DL is highly sensitive to any pa¬ 
rameter change, and is observed to have varied by a few 
millimeters between electrostatic probes and LIP mea¬ 
surements). Although LIP is efficient for negative ion 
fraction measurements in presence of an ion beam, this 
technique also has limitations. All necessary experimen¬ 
tal conditions were verified for all probe positions in the 
diffusion chamber (saturation of the collected current, 
colinearity of the probe tip with beam axis); however 
one should estimate the relative error to be around 15%. 
Another concern is the measurement of high negative ion 
fractions (typically n_/ne > 1). A typical temporal pro¬ 
file for the increase of the probe current (A/g) during 
photodetachment is shown in the insert of figure 1121' bl . 
The saturation value of the probe current (A/g used) 
was used for the determination of a according to equa¬ 
tion [5l However, according to Ref. [b^, a sharp peak 
may be observed at the beginning of the probe current 
pulse for large negative ion fraction. In this case the 
authors in Ref. |65l| recommend to use the peak value 
(A/g peak correction) in equation 0 Our signals show 
an overshoot at the beginning of the pulse and large high 
frequency oscillations are also observed. Since we are un¬ 
certain as to the origin of the high frequency oscillations, 
we ignored the peak correction. Taking this peak cor¬ 
rection into account would lead to a 50% increase in the 
negative ion fraction plotted in figure [T^b). Finally, one 
further limitation of the measurement is the uncertainty 
with the SF5 electron affinity, and the fact that the laser’s 
photons at 355 nm may not be able to detach SF^. An 
upper estimate of the SF^ fraction is 0.5 (Kono estimates 
this fraction to be a few hundredths in an ICP and 
was measured in another device to be 0.15 13). Thus, 
it appears that the actual negative ion fraction from the 
LIP technique lies between the measured value presented 
in figure dUb) and twice this value. 

Although no definitive conclusions on the absolute val¬ 
ues of the negative ion fraction could be drawn, the spa¬ 
tial evolution is clearly identified. The negative ion frac¬ 
tion experiences an abrupt increase at the DL position 
when entering the low plasma potential region. The DL 
separates a highly electronegative plasma in the diffusion 
chamber from an electropositive plasma in the source. 
Downstream of the DL, the negative ion fraction gradu¬ 


ally decreases from its highest value, at the DL position, 
to a very small value at the bottom of the chamber. This 
spatial evolution was captured in a recent model [b^ and 
seems to be a result of the plasma expansion. In the no- 
DL case, the negative ion fraction profile displays a shal¬ 
low maximum at the two chambers interface. However, 
either side of this maximum, the negative ion fraction 
decreases gradually. 

We believe that the value of the negative ion fraction 
measured slightly above the DL is actually zero. This is 
confirmed by the negative ion fraction measured by LIP 
in O2 presented in the following section. 

Though absolute values of the negative ion fraction 
published in Ref. [3l niight be overestimated, LIP con¬ 
firms that the DL is an internal sheath between a highly 
electronegative plasma (a > 1) in the low potential re¬ 
gion and a weakly electronegative plasma (a < 1) in the 
high potential region. Note that a weakly electronegative 
plasma has an electropositive pre-sheath and the nega¬ 
tive ion fraction could be as high as 1 in the center of the 
source. Unfortunately our experimental arrangement did 
not allow us to measure the negative ion fraction in the 
source. 


E. Parametric dependences 

The DL characteristics were investigated as a function 
of the external parameters. Figure [T5I al shows the power 
dependence of the potential obtained from the lEDF 
when operating at 1 mTorr, while maintaining the SFg 
concentration at 8±0.5% such that the DL position is 
at 25±1 cm. Figure [T5I bl shows the same analysis as a 
function of pressure in the range 0.5-b mTorr while the 
power coupled to the plasma is kept constant at bOO W. 
Variation of the coupled power does not influence the 
potential drop across the DL, which is consistently mea¬ 
sured to be 8±1 V from the lEDF analysis. The influence 
of pressure is somewhat surprising. DLs were observed 
for a broad range of pressure, from 0.5 mTorr to 10 mTorr 
- however above b mTorr, static DLs were only observed 
at very high power (above 2000 W) and propagating DLs 
observed at low power. As the pressure is increased, the 
particle balance imposes that the electron temperature 
decreases. lEDFs analysis shows a decrease of both the 
plasma potential and the beam potential while the poten¬ 
tial drop across the DL remains at 8±1 V when varying 
the pressure from 0.5 to b mTorr. An important observa¬ 
tion is the pressure dependence of the SFg concentration 
for DL formation: decreasing from 15% at 0.5 mTorr to 
5.5% at b mTorr. The independence of the electroneg¬ 
ative DL potential drop with pressure is distinct from 
the electropositive case, whose potential drop has been 
shown to be proportional to the electron temperature. 

As in the electropositive magnetized case, the influence 
of changing the boundary conditions has been investi¬ 
gated. We observed that for the same gas mixture, pres¬ 
sure and power, a change in boundary conditions does 
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not produce a significant change in the plasma potential 
profile. From geometry (a) to geometry (c) (presented on 
figure [7|), the plasma potential has the same spatial evo¬ 
lution, with a 2 V offset for geometry (c). The potential 
drop across the DL and the DL position was not modified 
by varying the boundary conditions. The formation of a 
DL in geometry (c) attached to the source tube forces 
the DL to be current-free, and we believe that any of the 
three geometries will lead to a current-free DL. 


F. The Ar/02 mixture 

DLs were also observed when operating in Ar /02 gas 
mixture. As for the Ar/SFg gas mixture, the DL forms 
when the O 2 concentration is above a critical threshold. 
Figure [M] shows the spatial evolution of the plasma po¬ 
tential for 68 % (dots) and 72% (open circles) O 2 concen¬ 
trations showing the DL formation at an operating pres¬ 
sure of 1 mTorr. For any O 2 concentration when diluted 
in argon, the plasma remains stationary, and the dis¬ 
charge does not enter the propagating DL regime, which 
is the main difference from the Ar/SFg mixture. The 
characteristics of the DL are similar to the Ar/SFg case: 
the DL separates (i) a high plasma potential region in the 
source and a low plasma potential region in the diffusion 
chamber, (ii) a high electron density plasma upstream 
and a low electron density plasma downstream, and (hi) 
a high electron temperature plasma upstream and a low 
electron temperature plasma downstream. LIP measure¬ 
ments of the negative ion fraction in Ar /02 gas mixtures 
is simplified n comparison to Ar/SFg gas mixtures: all 
negative ions created ( 0 “, O/" and O 3 ) can be pho- 
todetached by 532 nm photons. Thus, the measurements 
do not present a large uncertainty in the absolute values. 
The laser beam energy dependence of the LIP signals pro¬ 
vide a measurement of the photodetachment cross sec¬ 
tion and therefore of the dominant negative ion. For 
the present operating conditions (1 mTorr - 750 W), the 
dominant negative ion is 0 “ as is commonly observed in 
similar discharges Figure ITS] shows the axial pro¬ 

file of the negative ion fraction in the diffusion chamber 
for three O 2 concentrations ( 68 , 72 and 100%). Strong 
similarities with Ar/SFg gas mixtures are observed; the 
negative ion fraction profile is continuous with low ab¬ 
solute values for the non DL case, and for the DL case 
an abrupt drop is observed at the DL position. Absolute 
values of the negative ion fraction downstream of the 
DL are lower than in Ar/SFg gas mixtures due to lower 
O 2 attachment coefficients. These profiles show that the 
negative ion fraction is close to zero in the source region, 
which is consistent with the results in Ar/SFg mixtures. 
The negative ion fraction also experiences a decrease at 
the bottom of the diffusion chamber, which i s ty pical of 
stratification in electronegative plasmas [251 - 1^ . If we 
were to take into account the peak at the beginning of 
the pulse, the negative ion fraction would be increased 
by 20 % (see insert in figure fTSl). 


V. CONCLUSION 

An extended experimental study of static DLs in a de¬ 
vice composed of a source chamber attached to a larger 
expanding chamber has been presented. Positive ions cre¬ 
ated upstream are accelerated across the DL, thus leading 
to a positive ion beam in the downstream region. The 
positive ion beam could be utilized as a long term high 
power space electric propulsion system as was proposed 
by Charles Q- 

With an electropositive gas, the DL forms at a very 
low pressure if a highly diverging magnetic field forces a 
strong plasma expansion. The typical voltage drop across 
this electropositive DL is 20 V at 0.17 mTorr. Several 
boundary conditions were investigated at the end of the 
source tube and showed that the DLs are not necessarily 
current-free. They are strongly attached to the magnetic 
field lines and appeared to be very stable. In particular 
their position was not modified by the intrusive electrical 
diagnostics approaching from downstream, but crossing 
the DL with Langmuir probes would destroy it. 

DLs were also formed when using electronegative gas 
mixtures (Ar/SFe and Ar/ 02 ). In this case the magnetic 
expansion is not required and no restrictive conditions in 
pressure were observed. The DL formation appears to 
be related to the presence of negative ions and does not 
depend on specific atomic or molecular processes. The 
potential drop across the DL is independent of the elec¬ 
tron temperature and is about 8 V. Unlike the electropos¬ 
itive DL, the electronegative DL is a loose structure the 
shape of which is extremely sensitive to either the experi¬ 
mental parameters (gas mixture, pressure...) or intrusive 
diagnostics. Moreover, the electronegative DL becomes 
propagates at higher negative ion fractions. Finally, the 
electronegative DL could be traversed by the Langmuir 
probe without being destroyed which may be a conse¬ 
quence of its loose structure. 
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FIG. 2: Normalized RFEA IV curve (solid line), experimental 
lEDF (dashed line) and gaussian fit of the experimental lEDF 
(dotted line) showing evidence of an electropositive DL at 0.1 
mTorr, 90 G and 250 W. 


FIG. 1: Schematic of the experimental device and normal¬ 
ized axial magnetic field magnitude (solid line) and derivative 
(dots) used for investigation of the magnetized electropositive 
DL in section imi 
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a) Position (om) 
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b) Position (om) 


FIG. 3: (a) Local plasma potential (□) and beam potential 
(•) as a function of position, normalized lEDF are plotted in 
the insert and (b) local ion density (□) and beam ion density 
(•) as a function of position, and theoretical spatial evolution 
assuming a flat profile (dashed line) and a hemispherical pro¬ 
file (solid line) for the electropositive DL. 0.17 mTorr, 90 G, 
250 W. 



FIG. 4: (a) Schematic of the electropositive DL profiles cho¬ 
sen to fit the experimental data and (b) modification of the 
electronegative DL profile as a function of SFe concentration. 
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FIG. 5: Local plasma potential, beam potential and potential 
drop across the electropositive DL (a) as a function of the 
static magnetic field amplitude and (b) as a function of the 
pressure with the RFEA at z = 20 cm. 



FIG. 6: (a) Time resolved lEDF when the plasma is pulsed 
and (b) the associated beam and plasma potentials and am¬ 
plitudes as a function of time (zoom of the ignition phase) 
showing ignition of the electropositive DL. 
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FIG. 7: Boundary conditions terminating the plasma stud¬ 
ied and associated lEDF measured in the diffnsion chamber, 
downstream of the electropositive DL. 



c) z (cm) 

FIG. 8: Spatial variation of (a) the plasma potential, (b) the 
electron density and (c) the electron temperature in Ar/SFe 
mixtures: 4% SFe concentration (no DL - dots) and 8% con¬ 
centration (showing electronegative DL formation - open cir¬ 
cles). 600W, 1.5 mTorr. 
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FIG. 9: (a) The spatial variation of the eepf in the presence 
of the electronegative DL, and (b) the eepf upstream of the 
DL (z = 28 cm, black circles), downstream of the DL (z = 22 
cm, gray circles) and shifted upstream eepf (solid black line). 
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FIG. 10: An experimentally measured lEDF showing two 
peaks and the associated fit used for systematic analysis of the 
electronegative DL. Insert: lEDF measured when the RFEA 
is in the radial position. 
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FIG. 11: a) The local plasma potential (□) and the beam 
potential (•) as a function of position (lEDF are plotted in 
the insert) and (b) the local ion density (□) and the beam 
ion density (•) as a function of position, and the theoretical 
spatial evolution assuming a flat profile (dashed line) and a 
hemispherical profile (solid line) for the electronegative DL. 
600 W, 1.5 mTorr. 
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b) z 

FIG. 12: The spatial evolution of the negative ion fraction 
measured (a) by the two probes technique for two negative 
ion temperatnre estimates and (b) by aser induced photode¬ 
tachment (LIP) for two SFe concentrations. The insert in 
fignre (b) shows the increase of the electron cnrrent dne to 
LIP. 
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FIG. 13: The evolution of the plasma and beam potentials as 
a function of (a) power at 1.5 mTorr, and (b) pressure at 600 
W, with the RFEA at z = 10 cm. 
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FIG. 14: The spatial variation of the plasma potential in 
Ar/02 mixtures: 68% O 2 (no DL - dots), 72% O 2 (electroneg¬ 
ative DL case - open circles). 
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FIG. 15: The spatial evolution of the negative ion fraction 
from LIP for a 68% O 2 mixture (no-DL) and 72% O 2 and 
pnre O 2 mixtnres (electronegative DL). The insert shows the 
increase of the electron cnrrent due to laser induced photode¬ 
tachment. 





